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This paper deals with the analysis and design of the optimal robust controller 
for the fuzzy parametric uncertain system. An LTI system in which 
coefficients depends on parameters described by a fuzzy function is called as 
fuzzy parametric uncertain system. By optimal control design, we get control 
law and feedback gain matrix which can stabilize the system. The robust 
controller design is a difficult task so we go for the optimal control approach. 
The system can be converted into state space controllable canonical form with 
the a-cut property fuzzy. For optimal control design, we find control law and 
get the feedback gain matrix which can stabilize the system and optimizes the 
cost function. Stability analysis is done by using the Kharitonov theorem and 
Lyapunov-Popov method. The proposed method applied to a response of 
Continuous Stirred Tank Reactor (CSTR). 
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1. INTRODUCTION 

To design a proper controller for real time nonlinear system have many problems such as uncertainty, 
disturbance, unknown exact mathematical model etc. Uncertainty is either structure or Unstructured. 
Mostly the nonlinear system defined in terms of mathematical model doesn't have the exact parameter the 
calculated parameters are considered. So that type of model consists the parametric uncertainty. And because 
of that uncertain parameter, we have some information loss, which may be incomplete, unreliable. Also, the 
uncertainty affects the system response. So, we must design a controller such that it can deal with this type of 
uncertainty and remove the effects occurred because of that. For that, the fuzzy logic controller gives the best 
solution. In this paper, the design a controller which deal with parametric uncertainty with a fuzzy controller 
is given. That is the system with known mathematical model represent in terms of fuzzy membership function 
i.e. fuzzy parametric uncertain system. That means we consider some range of parameter variation. Then find 
the a-cut set for that fuzzy parametric uncertain system. From that we get the nominal system, then find the 
feedback gain matrix for that system which stabilize the system. We can design a control law for worst case 
condition i. e. for a = 0 when the greatest uncertainty Considered. The stability is checked by 
Kharitonov theorem. 

Robustness is the property of the system which deal with the variation in parameters in some range or 
bound [4]. In [5] this the method for fuzzy parametric uncertain system we discussed to design the robust 
controller with considering two examples. In first example, two uncertain parameters considered and in another 
example three uncertain parameters considered checked the response for that. The arithmetic calculation for a- 
cut set property of fuzzy are given in[1]. The stability of system is checked by Kharitonov theorem [3]. Because 
for this theorem we consider some range of parameter and in fuzzy parametric uncertain system uncertain 
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parameter considered in some range. In [2] the design of Popov-Lyapunov stability method is discussed. 
The design of paper is as follows, in section II represent the proposed methodology in that we discuss the a- 
cut for fuzzy parametric uncertain system. In section III we study the CSTR and apply the proposed method to 
for system and study the result. In section IV the stability of system by applying method is checked using 
kharitonov polynomial and Popov-Lypunov method [2]. In section V we end the paper. 


2. PROPOSED METHODOLOGY 
2.1. a-cut set 

Fuzzy parametric uncertain system is the co-ordinates of system depend on parameters represented as 
fuzzy membership function. A parametric uncertainty is present when system represent in mathematical model 
form, which is not exactly known. Here, the uncertain parameters can be defined in-terms of fuzzy function 
p, € p with membership function a = u(p;) € [0,1] as shown in Figure | The p is a universe of discourse 
for, ~;.The membership functions w(G,.) are having single mode and decrease to the interval last limit. Thus, 
the fuzzy uncertain parameter p, with a-cut is given by 


p,(a,) i [Dr (@,), De (a,)] (1) 


B@) 


pr Ca) Pr Pita) vet 


Figure 1. @ cut for fuzzy uncertain parameter p_i 


Where a;is the membership height for the, f, also p; (.) is a growing function and p; (. ) is a reducing 
function. Therefore, 


Pri(0) = PriPr (0) = PriPri (1) = Pi=Pri (2) 


The membership value a;can be taken as the confidence degree form that we get the nominal system. The value 
a, = 1 indicates the precise knowledge. 


Pri= ker (p;1) =Pri (3) 


Whereas a, = 0 represents maximum uncertainty 


Pri= supp (6,1) =[Pri. Pri] (4) 


2.2. Optimal controller design 
Consider a plant with uncertainty: 


Gn s+ + ...+G9 


#(s',a’,p") = (5) 


s'N4+Hy_ 451-14 ...4D19 


Where a’ , and p’ rt =0, 1, 2, ---,n-1, interpreted as fuzzy function. If a degree of confidence a € [0, 1] in 
the coefficients is given, then in (5) can be interpreted as an interval system as: 


Gt_4(a)si—1+4 ...4@79 (a) 
SIN+Din—1(a@) S14... + P79 (@) 


T(s',a'(@),p'(@)) = (6) 
Where 4,(a) and 4,(@), r=0, 1, 2, ---, n-1, interpreted as fuzzy function. The FPU system in (6) can be 


represent in state space controllable canonical form: 
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0 1 tee 0 0 
vel ga en fee fie 
—po(a) —fo(a) + —fip(a) 1 
y'=[-Go(@) -G,(a) + —G,-1(@)] (7) 
Where a € [0, 1]. The system in (7) represented as 
x = A'(p(a))x + Blu 
y’ = C'(G(a@))x (8) 


Where p’ is the parametric uncertainty vector in terms of a, and A, B, and C are the state space matrices. 
The uncertainty in state matrix (A’) represented by (8) is balanced with B’, if this uncertainty is in the bound 
of B’. Now a nominal value, Pyom € P(@) such that (A’ (Pnom), B’) is stable. Also, it is assumed that for any 


p' (a) € p there exists a 1xn matrix O( p’(@))so that, we can represent the uncertainty in State matrix (A’) as 
A'(p'(@)) = A' nom) = B'O( p'(@)) (9) 


Where @( p(a)) is bounded and Pygom € P(a@)is the nominal value of p(a) 


A' (p'(@)) = A'Pnom) + BOC P'(@)) (10) 
From (10) the (8) becomes, 

x! = A'Dnom(a@)x' + B'O(p'(a)x' + B'ul (11) 
is stable for all p(@) € p forall a € [0, 1]. Robust controller design for this problem is difficult. Hence optimal 
controller is designed for FPUS in (8). If the system is stable then only we can design LQR controller. The 
system in (9) is always controllable, hence it implies that pair A’ (p’(a)B’) is stable for each a € [0, 1] 
LOR design 


The problem can be solved using LQR controller which can be designed as: 
For the nominal system correlate to each a € [0, 1]. 


x= A(dnom)x' — Bu’ (12) 
A control law, u’=-k*x’ which reduces the cost functional 

[= fore +x!" x! + ul R'u')dt (13) 
Where F’ is an upper bound on the uncertainty O(p' (a))" 0(p'(a)), that is, for all B(a)ep 

0'(p'(a))"0(p'(a)) < F’ (14) 


The theorem in (11) can be used to design an optimal controller for FPUS in (8). When O(p' (a@))is enclosed 
the existence of the upper bounce on F’ is guaranteed. The cost function J can be written as 


J'.= Fe" Oe + ul R'u')dt (15) 


Where Q’= [F’+]] and R’= [1]. Thus, if (14) satisfied then an optimal controller designed. 
Error weighted matrix 

Design a controller for a=0 which consist of maximum ambiguity that controller stabilizes the all 
other system for different values of a. For « = 0, the parametric ambiguity is described by p, € [p; p;], for r 
=0, 1,2, ---,n—1. Calculate the error weighted matrix which satisfy the condition in (14) 
Consider the nominal value as, 
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P noni = [Po Pr wees Pn-1] 
Take the nominal value, 
P'nom = [Po Pi ----Pa-1] 


or any value in between [p; p;*] Then the nominal system of (8) is given by 


0 1 = 0 0 
P= |e 4 a 4 Peele 
—Po —P1 oo —Pn-1 1 
The uncertainty can be written as 

0 1 ss 0 01[ Po — Po 
: : ‘ 3 = |i||) Pi = Pa 

0 0 oo 1 ~ 10 ; 

Po —~Bo Pi -P1 * Pn-1—~ Pn-1 LV Lpn-1 — Pn-1 


Hence the uncertainty satisfies the matching condition. The uncertainty is given by 
0' = [po — Bo Pi — BiPn-1 — Pn-1] 

Which is bounded by 
OT O =F" 


Using F’ optimal controller is designed which reduces the cost function J. 


3. CASE STUDY CSTR 


A real-time experimental setup for highly nonlinear tank is constructed shown in Figure 2. DAC is 
used to interface CSTR with the Personal Computer (PC). The overall system consists of a tank, pump, 
Rotameter, RTD, an electro-pneumatic converter (I/P converter), a pneumatic control valve, an interfacing 
DAC module and a Personal Computer (PC). The differential pressure transmitter output is interfaced with the 
computer using DAC module in the RS-232 port of the PC. Figure 3 shows the block diagram of a CSTR tank 
interfaced with PC. The pneumatic control valve uses air as an input and adjusts the flow of the water pumped 
to the CSTR jacket from cold water tank. This flow maintains the temperature inside the tank at the 
desired value. The temperature of the liquid inside the tank is measured with the help of RTD and is transmitted 
in the form of (4-20) mA to the interfacing DAC module with the help of temperature transmitter to the 


Personal Computer (PC). 


Figure 2. Experimental setup 
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Figure 3. Block diagram of CSTR 


In Figure 3 the block diagram for CSTR gives in that we see the three regions first computer in that 
we design a controller which connect to the system through the data acquisition system, third is process flow. 
After calculating the control algorithm in the PC, required control signal in the form of current signal (4-20) 
mA is transmitted to the I/P converter, which passes the air signal to .operate by this signal to produce the 
required flow of water in and out of the tank. 


Calculation of TF: 

For calculating Transfer function of CSTR cooling process, the step response is taken into consideration. 
The transfer function is calculated by using process reaction curve. The process has very large dead time and 
is highly damped. Therefore, the step response can be fitted into a simple first-order model with dead-time. 


Ale 9's 


gG’® = 
Ust+1 


Where, A=Process gain, 0’ =Dead-time, t’=Time constant Therefore, the transfer function of the process is 
given by 


0.12e725 


oS) ae ee 


The Transfer Function of Valve is 


0.112 


G' = ~——— 
(= gsya 
By using Pade’s approximation the second order transfer function is calculated as 


—0.12s + 0.12 


ONO) 


The state space matrices are given as 


, _ [-1.33 —0.667 , _[0.25 
44 = [95 0 | Bo 
C’, = [-0.16 0.32] D',= [0] 
Simulink design: * 


3.1. Proposed Technique 
System transfer function 


0.0056 


G'sys(S) = 3575635 40417 
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Consider two uncertain parameters which represent in fuzzy no. as first a=tri (0.03 0.05 0.07) and 
second is b=tri (1.2 1.5 1.7) and c=tri (0.2 0.4 0.6). The acut for this is a= [(0.02a+0.03), (0.07-0.02a)] b= 


[(0.2a+1.2), (1.7-0.2a)], and c= [(0.2a+0.2), (0.6-0.2a)]. For different values of a the values of a, b, c are given 
in table 


Table 1. a—cut set 
a=0 a=0.25 a=0.5 a=0.75 a=1 
[ [0.035 0.065] [0.04.0.06] [0.045 0.055] [0.05 0.05] 
[1.21.7] [1251.75] [1.31.6] = [1.35155] [1.4.14] 
[ [0.25 0.55] [0.30.5] _ [0.350.45] [0.4 0.4] 


Feedback gain matrix: K’= [0.8954 0.9168] 
The closed loop characteristic Equation for a = 0: 


[1.1168, 1.5168] + [2.095, 2.59] s+s’=0 
Kharitonov polynomials are: 

K’;= 1.5168 + 2.595 +s° 

K’= 1.5168 + 2.0955 + s* 

K’3= 1.1168 + 2.0955 + s* 

K’4= 1.1168 + 2.59s +s” 


Step Response for different values of a € [0, 1] 


State response plot 


State response plot 


State 1 


o 5 10 15 “0 2 4 6 8 10 
Time (seconds) Time (seconds) 
Figure 4. Step Response of states for Figure 5. Step Response of states for 
different values of a different values of a 


In Figure 4 and Figure 5 shows the state response for different values of a. In that, the response for a=0 gives 
the stable response for other values a. In Figure 6 the bode plot for interval plant is shown. 
Bode plot for interval plant: 
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Figure 6. Bode plot 


Algorithms: 


Initialize population size, number of design 


variables and termination criterion 


Identify best and werst solutions in the 
population 


Modify the solutions based on best and worst solutions 
Kes =X prt yal Eves)-|Epradl re lX; vores) | O00] 


* Yes Is the solution corresponding to X 9 


better than that corresponding to y | ? 


y 


Accept and replace 
the previous solution 


Y 


Keep the previous 
solution 


No 


Y 


Report the optimum solution | 


Figure 7. Flowchart for JAYA algorithm 


Satisfy stop 
criterion? 


Select Genetic 
Operator 


Crossover Operator: 
Select two individuals and swap 
a section of gene between them 


Mutation Operator: 
Select one individual and 
mutate the genes in it 


Figure 8. Flowchart for GA algorithm 
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Initialize no. of students (population), no. of subjects (design 
variables), termination criterion 


Calculate the mean of each variable 


Identify the best solution (i.e. x,) 


Modify solution based on best solution 


xi(s) — x,(s)+ r&x,(s)—T,x¢s)) Teacher 
Phase 
Is solution x’ better than 
solution x,;? 
Keep the previous Replace the previous 
Select two solutions randomly x’ and x/ 
Is solution x’ better than 
solution x! 
*, r , , * ’ , ’ 
x (s) — x, (8) +r(xc,(s)—x,(s)) X,(S5)<— x,(s)+r(x,(s)—x,(s)) 
: : : Learner 
Phase 


Is solution x better than 


solution x’? 


+ 


Is the termination criteria 
satisfied? 


Report optimum 
solution 


Figure 9. Flowchart for TLBO algorithm 


Define the solution space, fitness 
function, and population size. 


Initialise position, velocity, 
pbest and gbest 


If fitness (x) better than 
fitness (gbest) 


Evaluate fitness function 


= : If fitness (x) better than 
Eger = ii (w — |D(e*)|)? dw fitness (pbest) 


No iteration = max. no. iteration 


or 
Fitness (gbest) is good enough 


Yes 


Solution is gbest 


Figure 10. Flowchart for PSO algorithm 
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f Initialize food sources Xj; } 


( Evaluate the solution quality for each X;; ) 


Initialize new food sources for employed bee 
and evaluate solution quality by these sources 


Initialize new food sources for 
onlooker bee and evaluate 
solution quality by these 


Compare solution quality 
by the new food sources 
and the initialized ones 


sources 


Determine food 
source Xj; with 
higher P, 


Are all onlooker 
bees distributed? 


Calculate position P, 
for food source X;; 


Yes 


Determine blocked 
food sources 


No 


Create new food sources to 
replace the blocked ones 


Evaluate total honey (best 
solution) for these food sources 


Store the position of the best food source 


Yes 


No 


Figure 11. Flowchart for ABC algorithm 


To avoid membership function as too redundant or too separate we use this algorithms for 
optimization. Fuzzy JAYA gives best result among all this algorithms namely GA, TLBO, ABC, PSO. 


20 T T T T T T T T T 
JAYA 
——GA 
15 —--—-TLBo LH 


Sg | ! 1 ! L ! 1 L L 
Oo 200 400 600 800 1000 1200 1400 1600 1800 2000 


Figure 12. comparative results by using different algorithms 


Table 2. Comparative Parameters with different algorithms 
Algorithms Rise Time Settling Time Overshoot Peak ISE 
JAYA 1.9947 4.5821e+03 1.0207e+04 18.0752 381.2 
GA 0.3289 5.052 1e+03 2.8753e+04 9.9626 415.7 
TLBO 1.0061 5.1323e+03 1.5342e+05 10.8502 475.7 
ABC 1.4707 5.582 1e+03 1.8207e+04 18.0752 475.7 
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PSO 1.4707 5.569 1le+03 1.8207e+04 10.5219 345.7 
3.2. Hardware Result 
For this, we have considered the error as 0-5 i.e. the difference between the setpoint and the actual 
value is 0-5 and control action is taken between 0-100. 


Time response plot 


., swenget point 
80k 4 — actual value 


Ss =-=control action 


Data 


“ 1 1 1 1 ! 1 1 
0 50 100 150 200 250 300 350 86400 
Time(sec) 


Figure 13. Hardware Result 


The figure. 9 shows the response for the real-time system. In this, we see that the set point is tracked 
by the actual temperature value. Also at the same time control action is also plotted. 


4. STABILITY ANALYSIS 
4.1. Kharitonov Theorem 
Consider a family n of real rational polynomials 


A(s) = Ays + Ags?+...... +An-1S" 1+A,s" (20) 
This polynomial family is called the family of uncertain or interval polynomial. 
xi< Ai< yi, where, i=0, 1, 2, ---, n. 


In accordance with Kharitonov theorem every polynomial in the family A(s) is Hurwitz if and only if the 
following four extreme polynomials are Hurwitz [3]. 


K', = Xo + x18 + yos? + y387 + 
K'p = Vo t+ y1S + X87 +.x387 + 
K's = Yo +xX18 + X28? + y3s32 + 
K', =Xo + VS + 2s? +.x387 + 


To apply this theorem, we need consider some range and for fuzzy parametric uncertain system we consider 
uncertain parameter in the form some range. 


4.2. Popov-Lyapunov Method 

The stability of fuzzy control system is checked by Popov-Lyapunov approach [2]. In order to do this, 
we transform the fuzzy system into Lure system with uncertainty. After that Lyapunov direct method is used 
to guarantee the stability than the robustness measurement which gives the bound on allowable uncertainty. 
The allowable bounds can help us to estimate the robust stability of fuzzy control system. The dynamic 
Equation of the system is 


a aay ee [al - (ile 
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re X14 
ee [1 0] [| 
From (21) we know theAg , Bo, Co. Now, first consider fuzzy control system, 
Where, r’ is reference input and wu’ is controller output. For this system. K,,. K, And. K,, are the gain values 


which will be 0.8954, 0.9168 and 1 respectively. And K, Kz , are the feedback gain values. 
The system can be representing in the form of, 


A’=A') —K',K' minb' oC’, , b' = —K',b'o [23] 
C’=-C,',% = A'x — b,po = —C'Tx the values of A’, b’ and c’ are, 
po Si a if i [OT ak 
4s=[o59 —-ssi96h? =h)¢=0 [24] 


The transfer function can be obtained as 
G’ = C(sI — A)1b [25] 


This is the Lur’e perturbed system. The system defined in (25) satisfies the following condition then the system 

is asymptotically stable 

1. The nonlinearity o always belongs to the sector [0, K,” ] whereK,” is a positive number. 

2. The system matrix A, is Hurwitz (G(s) is stable), and there exists a scalar 7 > 0 such that -1 / r 6=A;, where 
A;is an eigenvalue of A,, and 


1/K” + R'e[1 + jwr]G’ (jw) > 0,¥w € R’ [26] 


1 
3. Let = ar, 


positive-definite matrix, there exists a > 0, a vector q, P is symmetric positive-definite matrix and Wo, 
and 6 > 0 satisfying 


t+co,y=rC'bt+ — y= rc’ (27) where r is chosen such that y > 0. Where W is symmetric 


A'Tp' 4 P'A' a —q'q'" _ eW'P'b’ = v'Jy'd' 
=eW’ 
=eW'A'+6171' [28] 
The robust stability is found by satisfying theorems 1, 2 and 3 


1 


GS) = SF 1819s +259 


Time Senes Plot: 


o O2 0.4 O. 


6 08 1 12 
Time (seconds) 


Figure 14. Popov plot 
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In Figure 10 the Popov plot is shown. From this plot, we can find the slope of the line which is used for the 
further calculations. From the Popov plot for G Ga@), we get the value or r=0.5. 


VT oisos |and 
y = 5.0198 
Let eW’ = ea P| 


After solving the algebraic Riccati (28) get the P matrix as 


4.4455 0.5796 


B = levee 1.1618 


| and B=0.2591 [29] 


The P matrix is positive definite, so the system is stable. 


5. 


CONCLUSION 
The controller is designed for FPUS using the a-cut property of fuzzy set. Controller designed for a 


critical condition that is for maximum uncertainty interval and that will stabilize the other interval of 
uncertainty. This technique is applied to continuous stirred tank reactor and study the responses. Also, the 
proposed technique is applied with a fuzzy optimized membership function using Jaya algorithm. The stability 
of the system is checked by Kharitonov polynomials and Popov-Lyapunov stability theorem. 
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